function c = RHS(x,y)
c=[sum(y);sum(x.*y);sum(x.^2.*y)];